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We discuss on very general grounds possible lineshapes of composite particles with 
one unstable constituent. Expressions are derived in a coupled-channel formalism 
for constituents interacting in an S'-wave with no assumption made on the width of 
one of them. We show how easy-to-use formulae, already existing in the literature, 
emerge for narrow constituents and identify the parameter that controls the regime 
of their applicability. 



I. INTRODUCTION 



In recent years charm spectroscopy had a renaissance as a result of a huge amount of 
data that came from modern experiments. Many new states were found at 5-factories and 



Tevatron, with the bulk of them showing 
from quark models of the Cornell type 



properties in vast conflict with what was expected 
l| or variants thereof — for a recent review see 



Ref. [2|. The experimental situation is expected to evolve further with BES-III just having 
started data taking and with LHCb and PANDA on the horizon. 

If there are no nearby thresholds, the lineshapes of resonances may be presented in the 
standard Breit-Wigner form. However, the proximity of a (relevant) threshold distorts the 
lineshape strongly, as was stressed recently in studies of the X(3872), and then the lineshape 



of a resonance may contain important information on its nature 0, 4] - - under certain 
conditions, discussed below, this enables one to distinguish between genuine quarkonium 
states and states with a significant admixture of a hadronic molecule in the wave function. 



However, it was found in Refs. js, 5, 6] that the resonance lineshape might well get distorted 
significantly by the width of the constituents. In this paper we investigate on very general 
grounds how this distortion depends on the width of the constituents and under which 
circumstances simplified formulae for the line shapes can be used. 

Indeed, the list of new near-threshold states includes, together with the above-mentioned 

n 

X(3872) particle, located close to the DD* threshold, such entries as the F(4260) |7| (which 
is not far from the DD\ threshold, where D\ is a P-wave 1 + charmed meson) and the 
F(4660) state (close to the ^7o(980) threshold). While the width of the D* is tiny (less 
than 100 keV), constituent widths, generally, can be as large as 50 MeV for / (980) or even 
larger. In this regard, we address the following questions: 

• What are the most general expressions for the lineshapes with no assumption made 
on the width of the constituents? 

• What is the parameter that controls the deviations from the simplified expressions of 



lU fs. 

The formalism employed is a straightforward generalisation of that used in Ref . {9} , where 
only the case of stable constituents was discussed. We believe that our insights will prove 
useful for the analysis of future experiments. 



II. COLLECTION OF MAIN RESULTS 



In order to simplify the reading of the paper and to allow those interested in the final 
result only to skip all technical details, we collect in this chapter the main results and 
conclusions — the derivation of those is then presented in the subsequent chapters. Indeed, 
while the derivation of the main results turns out to be rather technical, the final conclusions 
are not. 

For simplicity we denote the decaying particle as X. Although the results derived are 
more general and apply not only to the X(3872) charmonium, we shall use the latter as a 
paradigmatic example for the application of the suggested approach. We use the following 
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shorthand notations: in the intermediate state there are particles a (stable) and b (unstable, 
with the width r R ); in the final state we have the same particle a, while the particle b 
decays into particles c and d. Thus, in these notations, the process under consideration is 
X — > ab — > a[cd] (the corresponding chain for the X(3872) obviously looks like X(3872) — > 
DD* — > D[Dtt]). The masses of the particles are labelled accordingly (see Fig. [T]). Then, as 
outlined below, in most relevant cases the differential rate close to the ab threshold, with a 
relative 5*-wave between a and b, may be written as [see Esq. fl69|) below] 

dBi(a[cd]) _ ^\ gabKs{E) 

dE ~ 2tt|£_£ x + l g ^[ KeS (E)+ik eS (E)] + ir \ 2: 
where the constant B absorbs all details of the X production, g ab denotes the effective 
coupling constant of the X particle to ab, and n e s(E) and k e g(E) parameterise the real 
and the imaginary part of the X self-energy. The energy E is counted from the three-body 
threshold, so that the full mass is M = m a + m c + + E. Similarly, the parameter Ex is 
related to the physical mass of the X particle as Mx = vn a + m c + + Ex- In addition 
Eq. (CQ) contains Jo to mimic the presence of extra, inelastic channels (such as 7t + tt~J/i/j, 
7T + 7r~7r° J/i/j, etc. for the X(3872)). The corresponding inelastic rate then reads [see Eq. ( 1701) 
below] 

^Br(inel) = fi J_ ^ (2) 

dE " 2 *\E-E x + lg ab iK cS (E)+tk eS (E))] + tr \ 2 ' 
For a vanishing width of the constituent, J 1 ^ = 0, Eq. (JTJ) reduces to the well-known Flatte- 
type expression with 



k eS (E) = J2fj, p (E - E R ) Q(E - E R ), (3) 



k cS (E) = -y/2^ p (E R - E) &(E R - E) + ^2^ P (E R - E x ) B(E R - E x ) (4) 

for the relative momentum of the ab system and its analytic continuation below threshold, 
respectively. Here E R denotes the energy of the unstable constituent, also measured relative 
to the three-body threshold, via m h = m c + m d + E R , and \i v denotes the reduced mass of 
the ab system. Note that, deviating from the standard expression, K e fr was subtracted at 
E = Ex in order to comply with the definition of Ex as a zero of the real part of the inverse 
scattering matrix. One should keep in mind that all energies are defined with respect to 
the three-body threshold. However, in the limit r R = this threshold becomes irrelevant 
and the only relevant threshold is the ab two-body threshold located at E = E R — see the 
discussion below Eq. fl32|) . 
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Introducing a nonvanishing width r R implies a certain smearing of both K e g(E) and 
k e g(E), so that, in the presence of a finite width ^ 0, more complicated expressions 
for K e s(E) and k e s(E) need to be used, and we derive such expressions below [as given in 
Eqs. (14T|) - (I50]) ]. One might expect that for an isolated, narrow constituent, distant from any 
threshold, the situation simplifies, and this is indeed what we observe. We therefore define 
the dimensionless ratio 

with Er being the energy of the unstable constituent measured with respect to the three- 
body threshold. In this work we only study a resonance coupling to one channel — thus 
the constituent is isolated by construction. The parameter A then captures the other two 
conditions, since for A 1 the width of the constituent is much smaller than the distance 
to the nearest threshold. It will be demonstrated below that for this situation much simpler 
formulae emerge — our exact expressions reduce to those presented in Ref. \^ long time ago 
and, more recently, in Ref. [3|, namely 



k<§ (E) = ^jrj a/ (e - E R y + T 2 /4 + E-E R , (6) 



= -VJ^yViE - E R y + rV4 - e + e 



R 



+^yV(E x - E R f + rV4 - e x + e r , (7) 



where, analogously to above, a subtraction was introduced for k^. In Ref. jlj] the width 
was taken as constant, r = On the other hand, in Ref. [jj the same expressions were 
given, but with a "running" width, 

/ e \ ,+u - 

nE) = r R (-) , (8, 

where / is the orbital angular momentum in the subsystem cd. We shall demonstrate below 
that for both S- and P-wave 1 resonances as constituent, in the case of A 1, Eqs. ([6]) and 
(J7|), when used in Eq. (OQ), describe the lineshapes quite accurately. In this limit the width 
of the resonance can safely be considered constant and thus there is no improvement to use 
an energy- dependent width from Ref. instead. However, we shall also show that already 
for values as small as A ~ 0.2 significant deviations become observable between the full 
expressions and the approximate ones. Especially, for the approximate expressions, and this 
holds for both Eqs. ([6]) and ([?[) with a constant as well as with an energy-dependent width, 
the amplitude shows a wrong behaviour near the three-body threshold, which is part of the 
signal region for values of A ~ 0.1 and larger. The results for the lineshapes are illustrated 
in Figs. EHS 

We also discuss the effect of a possible interference in the final state. For example, since 
the X(3872) decays in equal parts via D°D°* -> D [D°n ] or D°D°* -> D [D n% the final 
states of both decay chains interfere. We argue that, to avoid a violation of unitarity, one is 
to include interference effects from the very beginning, by solving the full coupled-channel 
problem. However, for the constituent width small enough, a simple modification of Eq. ([T]) 
is possible which accounts for the interference effects, namely 

dBi{a[cd}) _ B l g a b(k c $(E) + k- mt (E)) 

dE ~ 2vr (E-E x + 9 -fK cS (E))i + | (P + g ab (k eS (E) + k int (E))) 2 ' 

For a motivation for this expression see text near Eq. (179|) . Notice that the interference 
contribution feint C^)) defined in Eq. ( J77|) . is taken into account both in the numerator and in 
the denominator of the differential rate as demanded by unitarity. We call this "improved" 
interference — to be confronted with the "naive" interference discussed in the literature, 



1 Notice that the P-wave case is relevant for the X(3872), assumed to be a D*D bound system, with the 
D* decaying into ttD in a P-wave. 
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with the fei nt (E) term included in the numerator only [compare Eqs. (179|) and (!7H|) and see 
Fig. HO]. 

This finalises our collection of main results. Before discussing them in further detail we 
now first turn to their derivation. 



III. GENERALITIES 



We now start with the detailed derivation of the results presented above. For simplicity 
we assume that the dynamics in all subsystems can be treated nonrelativistically. Our Fock 
space consists of three states: a compact seed, labelled as \X Q ), the ab pair, and the a[cd] 
state. Thus the wave function of the X can be written as 



( C\X ) \ 



(10) 



X(p)\ab) 
\<p(p,q)\a[cd\) J 

Here p and q are the centre-of-mass momenta in the ab and cd subsystems, respectively 
see Fig. [TJ 

Wave function ( TTTJj) obeys a Schrodinger-like equation: 



with the Hamiltonian 



If we label the channels as 



H\V) = Af|tf), 



H = H + V. 



(11) 



(12) 



1 = \X ), 2 = \ab), 3=\a[cd}), 



(13) 



we may write 



(o) , P 



p 2 q 2 



2// 



ab 



H 03 = m a + m c + m d + h — 



with the reduced masses: 



m a ml 



(0) 



m a (m c + m d ) 



m c m d 



m a + m c + rrid ^ q m c + m d 



m a + ml 

The masses and entering expressions (|14p and ffl5|) are bare masses. 



(14) 



(15) 
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Transitions allowed are \X ) -h- \ab), mediated by the transition potential V12 = V%\ = 
f(p), and \ab) <H- \a[cd]) (particle a being a spectator), mediated by the transition potential 
V23 = V32 = g(q)- Notice that, in addition to the transition potentials g(p) and f(p) 



one might allow for a diagonal interaction in the ab channel. It was shown in Refs. 10, 



111 ], however, that nonperturbative diagonal interactions in the mesonic channels can in 



general be removed within effective theories employing a unitary transformation — this 
is not possible only in the unusual situation of a hadronic molecule and a compact state, 
both weakly coupled, being present near threshold simultaneously 11]. In this way all 
nonperturbative meson-meson interactions get cast into the effective coupling constants 
and the parameter M^°\ As we shall see, the latter quantity does not play a role in the 
whole analysis, while the former contain important physics information — see Sec. |Vj 

We shall focus on resonances near threshold. This allows us to use nonrelativistic kine- 
matics and to identify /i a & = The latter implies that, in the expression for the reduced 
masses the bare mass is to be replaced by the physical mass ra& (to be defined below). 

Thus 



( M(°) f(p) ^ 



H = Hr ) + V 



(16) 



f(p) m a + mf ] + T ab (p) g(q) 
y g(q) m a + m c + m d + T ab (p) + T cd (q) J 

where two convenient abbreviations were introduced: 

T ab(p) = — , T cd (q) = — , 

for the kinetic energies of the ab and the cd subsystem, respectively. 

Now one can define the fully off-shell t-matrix as a solution of the matrix Lippmann- 
Schwinger equation, written symbolically as 

Uk = Vik — VijGji tik, (17) 

where the subscript labels channels and summation over repeated indices is understood. The 
free Green's function is 

g^ = (H -M)-\ g^ = 5 lk gl 0) . (is) 

For convenience we express all energies with respect to the thresholds, thus writing 
M = m a + m c + m d + E, = m a + m c + m d + E Q , mj; = m c + m d + E^°\ (19) 
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for the total energy E, the bare energy of the X, E 0l and the bare energy E^ of the particle 
b, respectively. 

The solution of Eq. (fTTl) takes the form 

G^(E)(E-E ) 

E-E, + Gf{E\ 

f( p)D 2 (p,E) E -E 
A{p,E) E-E + gP(E) 



ME) = Z XK :^^ (20) 

t 21 (p, E ) = t 12 ( P ,E) = fip ^zr„ :°'; ( o)^^ m 



A(p,E) 

D 2 (p,E) f( P )f(p') D 2 W,E) 
A(p,E) E -E Q + Gf{E) HV,E) ' 

ls g(q')D 2 (p',E) 



t 23 (p,p',q',E) = t 32 (p',q',p,E) = -5(p - p' , ^ y) 
g(q')f(p)f(p')D 2 (p,E) 1 



A(p,E)A(p',E) E-E Q + G ( °\Ey 



(24) 



+ ( n „ ' „> w\ „f\ 9(q)g(q') , g(q)g(q')f(p)f(p') i 

t 33{P , q , P , q ,E) = *(j>-j>)^^+ A(P , WB) E _ Eo + G f {E y < 25 > 

where the bare self-energy parts 

^-h ^ + T M -?-^m-« ' (27) 

drive the propagation of particle b and the X, respectively, and 

D 2 (p,E) = T ab (p) + Ej ) 0) -E, (28) 

A(p,E) = G i b °\p,E)-D 2 (p,E). (29) 

Let us now introduce physical masses m b and Mx and the corresponding energies E R 
and E x as 

m b = m c + m d + E R , M x = m a + m c + m d + E X - (30) 
We focus on the case of particle b being a resonance, thus having E R > 0. 
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The renormalised self-energies are then written as, employing subtractions at the physical 
masses (or energies E R and Ex), 

G b (p, E) = Gf (p, E) - Re (cf (0, E R )) , (31) 
G X (E) = [ d 3 p- f2{p) 



E R + T ab (p) -E- G b (p, E) - tO 

~ Re (/ £p E R + T ab { P )- f Ex P) -G b ^E x )- l ^) ■ (32) 

In such a way, the physical energies E R and Ex are defined as zeros of the real parts of the 
corresponding dressed propagators. 



Note that the self-energy of the X has a very peculiar analytic structure (see Ref. [12| and 
references therein): since G b has a branch point singularity at the opening of the cd threshold, 
the X self-energy shows two branch points inside the complex plane on the unphysical sheet 
- for the unstable resonance the physical sheet is the one where all three particles a, c, 
and d are physical; the two sheets meet at the three-particle threshold. To understand the 
role of these singularities better, it is instructive to observe the movement of those branch 
points as the width of the resonance decreases: then they approach the physical axis while 
at the same time the coupling to the physical sheet gets weaker and weaker. In the limit of 
a vanishing width of the resonance, what used to be the physical sheet decouples completely 
converting what was the unphysical sheet to the physical sheet; at the same time the branch 
points merge on the physical axis to produce the branch point that starts the unitarity cut 
for the ab channel. 

Equation (iTTj) for the wave function possesses three types of solutions (which can be 
defined by the corresponding asymptotic behaviour): (i) bound state of all three outgoing 
particles (a, c, and d), (ii) bound state in the subsystem cd and a free motion between 
this bound state and the particle a, and (iii) a free motion of all three particles, a, c, and 
d. Below we assume that neither a three-body bound state nor a bound state in the cd 
subsystem exists, and concentrate on the solution (iii), which is the situation most relevant 
for the X(3872) being located close to the D*D threshold with the D* as a resonance in the 
Dn channel. 

All relevant matrix elements can be found from the multichannel t-matrix discussed above. 
However, there are no prospects for experimental studies of, say, DD* elastic scattering, or 
of the inelastic process DD* — > DDtt. All information on the resonance X is therefore 
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extracted from its production from an external source. For example, from B — > KDDtt in 
case of the X(3872) particle. Then, if the said source can be considered as pointlike, the 
production amplitude may be written as 2 

M = constxg[ 0) (E)t l3 (p,q,E) = const x /(pM«) 



[E-E x + G X {E)\ [E-E R - T ab (p) + G b (p, E)} ' 

(33) 

where the details of the production dynamics are absorbed into the overall constant multi- 
plier and the momenta p and q satisfy the total energy conservation condition, 

E - T ab {p) - T cd {q) = 0. (34) 

For future convenience it is instructive to rederive Eq. fl33l) directly from the solution of 
Schrodinger equation (|TT|) . Indeed, for the case of a free motion of all three particles in the 
final state we may write 

^f-W, (35) 
where is the wave function of the free motion, 

^°\p,q\p,q) =5 i3 5(p-p)5(q-q), (36) 

in line with the assumed asymptotic. 

The solution of Schrodinger equation ( {Til is 

C(6 d ) = f(P)9(<i) (m 

{Pq) [E-E x + G x (E)][E-E R -T ab (p) + G b (p,E)Y 1 ) 

, i-x C(p, q)f(p)+g(q)5(p-p) 
XiPlP) = E — E R — T ab (p) + Gb(p, E) ' (38) 

t \~ ~\ u ~\xt ~\ i 9(q)x(p\p) i qc\\ 

V (j>, 9 \p, „) = Hp - p)% -i) + E _ TM _ Tcd(q) ■ (») 

Notice that the wave functions are expressed in terms of physical energies E x and Er. 



2 More correctly this is the expression for the production through the quark-antiquark component of the 
X wave function (see the discussion in Ref. 



13|) or if there is only a hadronic molecule present. It might be 



distorted significantly, if a weakly coupled molecular state and a quark state are present simultaneously in 



the near-threshold region — see Ref. 
this case further here. 



11] . However, since this is a quite exotic scenario we do not consider 
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It turns out then that, in practical applications, the quantity of an immediate relevance 
is the so-called spectral density, which can be built as 



W(E) 



d 6 pd A q \C(p, q)\ 2 5 (E - T ab (p) - T cd (q)) 



(40) 



and which measures the probability to find the bare \X Q ) state in the |a[cd]) continuum 
see Sec. [V] An explicit calculation gives 



W(E) 



2ni [E-E x + G* X (E) E — E x + G X {E) 
Comparing Eqs. (1331) and (1371) one can find for the differential rate 
rfBr (a [cd]) 



dE 



oc J d 3 pd 3 q\C(p, q)\ 2 5(E-T ab (p)-T cd (q)) = W(E). 



(41) 



(42) 



IV. SELF-ENERGY CALCULATION 



For energies E small enough one can neglect the form factors — then the resulting shapes 
show universal, mo del- independent patterns. In this limit we may write 



G b (p,E) = -g cd 



2fJL q [E- 



p- 



1+1/2 



e e- 



pr 



+ 2(-!)^ 



-2/i, E - 



2/-tp J \ 2fip 

e ( —E 



v 



2/x ? 



(43) 



where / is orbital angular momentum in the cd subsystem, and the coupling g cd is related to 
the physical resonance width .Tr via 



r R = gcd (^ q E R y +i/2 = 9i e 1 ; i/2 



J R 



(44) 



where, in order to simplify the notations, the shorthand gi = g cd (2ij, q ) l+l l 2 was introduced. 
Notice that, once we always stay in the vicinity of the resonance, we retain in Eq. ( 14*3]) only 
the part which changes appreciably in the region of interest, while a nearly constant (infinite) 
contribution is omitted in Eq. ([4*3]) as it can always be got rid of through an appropriate 
redefinition of the physical quantity E R . We should also stress that the case of / ^ is more 
subtle, as the form factor f(q) is proportional to q l , so that the renormalisation procedure 
requires more than one subtraction. However, as mentioned here we are primarily interested 
in the leading nonanalytical near-threshold contribution from the [cd] self-energy loop, which 
is given by expression (j43l) . 
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Thus, the two parameters we treat as input are the energy E R and the width r R . An 
important dimensionless parameter which governs the dynamics is 



A 



2E 



R 



(45) 
and shall 



We assume A < 1 — only then one can speak of the resonance as a constituent 3 
be primarily interested in the case A<1. 

We stick to zero orbital angular momentum in the ab subsystem. It is then convenient 
to express Gx{E) in terms of the coupling g ab as 



1 



1 



Re (G X (E)) = -g ab K eS (E), Im (G X (E)) = -g ab k cS (E), g ab = 8ttV p / 2 (0) 



(46) 



For a resonance in the cd subsystem one finds by an explicit evaluation of Gx{E), as 
given in Eq. (132|) . for the imaginary part 



kes(E) 



9i 



fdp(E - £- p ) (2l+1)/2 



V A \2l+\ 
2/V 



(47) 



^vh (E R -E + ^-y + ^(l 

where I denotes the angular momentum of the cd system in its rest frame and gi was defined 
in Eq. (jUJ). The real part is 

n eS {E) = K\(E) + k 2 (E) - k x {E x ) - k 2 (E x ), (48) 

where the first term takes the form 



K X {E) 

and the second term is 

K2(E) = - 



TT/^p JO 



9l 



p 2 dp- 



(^-£ + £) 2 + ¥(£-S) 2Z+l 



(49) 



p 2 dp- 



• 2^ p 



f_ _ E y2i+i)/2 



21+1 



(50) 



2lTfIp J ^ (E R -E + £- p y + £(E-£] 
Notice that K2(E) is suppressed as compared to K\(E) by a factor gi oc A, which is small by 
assumption. We neglect ^{E) in what follows throughout the paper. 
For example, for / = a straightforward evaluation gives 

/ \ 

g y/2fi p 



k cS (E) 



E R +\J(E-E R y + 



\ 



^E R (E R -E + ^ + x/(E - E R f f^pj 



(51) 



5 For a discussion of the scattering amplitude for general values of A see Ref. 



13 




2 4 6 8 10 12 14 

E[MeV] -80, 



FIG. 2: The exact results for k e g(E) and K e g(E) for the case / = versus the simplified formulae ([6]) 
and ([7]) with an energy dependent width as proposed in Ref. [3] for Er = 7 MeV and -Tr = 0.1 MeV 
(A = 1/140). The left plot is for k e g(E), the right plot is for K e s(E). The k% term is neglected. 



and 

k&{E)kk x (E)-k x (E x \ (52) 

with 



1 - £-^ + t-V(^-£) 2 + ¥ _ 

Ki(E) = const + yjHp — (53) 

'e r -E-§ + J(E R -E)* + & 



Notice that the (infinite) constant is cancelled in K e g, as follows from Eq. (1521) . 

For / > the integrand in K e ft shows, in addition to the two physical poles, spurious 
poles in the region p ^> ^2/i p E. Their appearance is an artifact of the formalism used 
[in particular, a consequence of neglecting the form factor in expression ( 13T|) for the self- 
energy of particle b]. However, since these poles are well outside the region of validity of the 
formalism discussed, in what follows we neglect their contribution (using the principal value 
prescription). Then analytical expressions for the k e g and k e r can be found, in analogy with 
the case of Z = 0. However, such expressions are bulky, so we do not give them here and 
resort to numerical calculations instead. 

To demonstrate quantitatively the comparison of Eqs. (117]) . (118]) versus Eqs. (JSJ), (JTJ) - 
for the sake of concreteness in all plots we use the energy-dependent width as proposed 
in Ref. — we choose two possible parameter sets. As the first example we choose the 
parameters as they apply to the case of the X(3872): for the width of particle h we take 
an empirical value of order of the D* state width, = 0.1 MeV. In addition, the mass 
of the D*° is located at about Er = 7 MeV above the D°tc° threshold. This corresponds 
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FIG. 4: The same as in FigEJ but for r R = 3 MeV (A = 0.2). 



to A = 1/140. The resulting energy dependencies of both k e g(E) and K e s(E) are given in 
Figs. [2] and [3] for S- (I = 0) and P- waves (I — 1), respectively. Note that the latter case is 
relevant for the X(3872), with the D* decaying into in a P-wave. Since we subtracted 
both expressions for K e s(E) at the same point E = Ex, the two curves agree exactly at the 
subtraction point, which is chosen here as Ex = Er — 0.5 MeV. Clearly, for the parameters 
relevant for the X(3872) the results agree perfectly well. However, already for a width of 
only about 1 MeV, deviations between the formulae start to become visible. In Figs. [Hand 
Owe compare the two formalisms for I# = 3 MeV (or A ~ 1/5) with all other parameters 
unchanged. Now the difference, especially in k e s, is sizable, mainly as a result of the very 
different threshold behaviour. Indeed, at the three-body threshold, Eq. (J171) behaves as 

k ^ E K%>J¥MT x E ' +2 • (54) 

in line with the three-body phase space and the centrifugal barrier. 

In the meantime, k^(E), as given by Eq. ([6]), shows a wrong threshold behaviour. If a 
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FIG. 5: The same as in Fig. [H but for I = 1. 



constant width is used for the resonance (this case is not shown in the figures), there remains 
a finite imaginary part even at the threshold, while for the energy- dependent width ([8]) of 
Ref. {2I one finds 

^ E ?o%fF(0«* ,2 ' +,,/2 < 55 > 

Thus, instead of the required behaviour according to the three-body phase space, the ap- 
proximate result scales as the two-body phase space. 

One concludes therefore that, if the resonance is well separated from the three-body 
threshold, the simple formulae ((6]) and (I7j) suggested in Ref. {5] may be used instead of 
Eqs. ( HTl) and ( HHl) . Indeed, it is easy to verify that, in the limit 

A«l, {Er ~ El « 1, (56) 

our formulae reduce naturally to Eqs. (j6j), ((7j) [it is an easy task to check this reduction 
analytically for the S"-wave case — see Eqs. (l51~j) and (|52|) ]. Notice that in this case the 
width of particle b can safely be taken as constant in Eqs. ([6]) and ([7]) and there is no 
improvement to use an energy- dependent width from Ref. 3| instead. This is definitely so 
for the parameters relevant for the X(3872). However, when investigating particles with 
broader constituents, one has to resort to the full formulae ( 14 7p and (14"9|) . This applies, in 
particular, in studies of the V(4260). The V(4260) resides at the D\D threshold; there are 
two Di states, a narrow one with a mass of 2420 MeV and a width of 25 MeV, and a wide 
one with a mass of 2430 MeV and a width of about 400 MeV, both decaying into D*ir 15 ]. 
In the former case one has A ~ 0.05, while in the latter case A ~ 0.8. 
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V. MODIFICATION OF THE WEINBERG FORMULAE 

The Weinberg formalism establishes a relation between low-energy observables, such as 
the effective-range parameters jlo| or Flatte parameters Q], and the nature of a resonance. In 
particular, it allows one to quantify the admixture of a genuine compact Xq component and 
a molecular component in the physical state X in a model-independent way. The formalism 
was designed to deal with two-body final states, and in this chapter we address the question 
as to how the formulae are modified by the presence of an unstable constituent. 

The central quantity is the spectral density which, in the near-threshold region, takes the 
form [see Eqs. (jUJ) and ( 1461) ] 

1 1 



W(E) = ~~~ 



(57) 



2ixi [E — E x + \g a b^cs — ^9abKs E — E x + \g a b^es + ^9abKs. 
which is nonzero above the lowest relevant threshold (in our case for E > 0). In the absence 
of three-body bound states it is normalised to unity, 



oo 



W{E)dE = 1. (58) 



o 



In the vicinity of the threshold at E — Er the integral 

rE R +8 

/ W{E)dE (59) 

JE R -5 

measures the admixture of the bare state for the near-threshold resonance. Indeed, if this 
integral appears of order unity, then one deals with a genuine quarkonium while a small 
integral value indicates a dynamical origin of the resonance. Therefore, fitting the production 
rate (j4"2|) with the formula (157|) . one can estimate the admixture of the bare X state in the 
wave function of the X-resonance in the near-threshold region defined by the parameter 5 
9j. Most naturally, the scale for S is set by the visible width of the resonance, while other 
scales, like the distance from the next threshold [D + D*~ in the case of the X(3872) particle] 
could also be relevant. 

In the narrow-width limit (1561) one can proceed further. Strictly speaking, with the 
finite width of particle b taken into account (and in the absence of three-body or two-body 
bound states), the physical space is spanned by the genuine three-body continuum, and the 
notion of the ab on-shell scattering amplitude ceases to exist. Nevertheless, for small enough 
energies and in the limit ( 1561) . the t-matrix element £22 takes the form 

t 22 ( E ) = — 9 M : , (60) 

8ll 2 fl p E - E X + ^abhs + 29abK e g 
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which can be viewed, for r R — y 0, as a Flatte-type representation for the t-matrix in the ab 
channel. 

With expression ( 1601) one can study the structure of singularities in the complex plane. 
If one of the particles is unstable, the Riemann surface does not show a two-body but a 
three-body cut (if the unstable particle predominantly decays via a two-particle decay, as 
is assumed throughout this paper). As shown above, in the narrow- width limit, one can 
replace the expressions for the k e s(E) and K e s(E) by the simple ones, as given by Eqs. (jUJ) 
and ([7j) with a constant width. Then the t-matrix poles are 



^gabKx - ~glbV P ± J Y^atb^l - \slbVp (e x -e r + ^g ab K X J - 



(61) 

with Kx = k^(Ex)- In the limit of a stable particle b, = 0, and for Ex < Ep>, there is 
a bound state in the ab system, which corresponds to the upper sign above, 

E 1 ^E B = E X , (62) 

and it is straightforward to verify that this pole is situated at the first sheet of the bona fide 
Riemann surface of the genuine two-body case. For r B > the system couples to the sheet 
with a, b, and c on-shell, and what used to be the physical sheet for stable constituents, is 
now the unphysical one. However, the pole E\ remains to be the one most relevant for this 
discussion. In its vicinity, £22 takes the form 

h2(E) -^L, (63) 



where 



gls 



gab 



g a b 



id. 



-2fi p (E B — e r + ~r A 



R 



8n2flp y/-2fi p (E B -E R + ir R ) + \ii p g ah 
In the limit ( l56l) and for E < Er one has k e s — > 0. Then 



(64) 



nE)^ , f EB J = 2S(E-E B ). (65) 

1 - ^gabQE^-^^pKti - tiR)\ E =E B 

The factor Z on the right-hand side of Eq. ( )65|) is nothing but the probability to find the 
bare state in the wave function of the bound state. 



18 



Equation ( 165]) allows one to express the coupling g a j, in terms of Z and Eb- Then the 
vertex g eS takes the famous Weinberg form [lo| 



As it was argued previously [10|, the given expression is correct only up to terms of order 
0{^/2fx p (E R - E B )/P), where denotes the range of forces. Our analysis reveals that 
in case of unstable constituents there are, in addition to the finite-range corrections, also 
corrections of order A. However, for A 1 the Weinberg analysis should still be applicable. 

VI. LINESHAPES 

We are now in a position to discuss the lineshapes of composite particles with unsta- 
ble constituents. For the plots shown below we consider the X(3872), as a paradigmatic 
example, which allows us to quantify also the effect of the inelastic channels in line with 
Refs. 0, 4|. 

We start with the expression for the spectral density, similar to the one given by Eq. (1571) . 



W(E) = — 
v ; 2ni 



1 1 



[E - E x + G* X (E) E-E x + G X {E)\ 



(67) 



with 

1 i 

G X {E) = -g a bK cS + - (gabKs + A) . (68) 

In this expression we introduced, in addition to what we had before, the quantity i~o to 
mimic inelastic channels (for example, the channel mrJ/ip in case of the X(3872)). 

Then the differential branching fractions are given by (as mentioned before, we absorb 
all details of the X production into the constant £>): 

dBi(a[cd]) = B l g ab k eS (E) 

dE ~ 2n(E-E x + afK cK {E)f + \{r G + g ah k c ^E)) 2) 

for the a[cd] channel (DDir channel for the X(3872)), and 

rfBr(inel) = I Jo (7Q) 

dE " 2-n(E-E x + 2fK e z{E)y + \{r + g ab k cS {E)f' 

for the inelastic channel. 

To illustrate the effect of the finite width of one of the constituents we investigate the 
above-mentioned differential branching ratios for both the S- and P-wave case and for both 
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FIG. 6: The a[cd] production rate [DDtt rate for the X(3872)] (the first plot) and the inelastic 

rate [ttttJ/^ for the X(3872)] (the second plot) for the S-wave resonance and for the X being a 

bound state. Here E R = 7 MeV and r R = (solid line), r R = 100 keV (dashed line), r R = 500 

keV (dotted line), and r R = 1 MeV (dash-dotted line). Notice that, since we are interested in the 

lineshapes, the above rates are plotted for B = 1. 
0.1 1 



;0.0'5 



-a 



a/- 



E[MeV] 



10 




FIG. 7: The same as in Fig. [61 but for the X being a virtual state. 



bound and virtual X particles. However, before we proceed, let us remind the reader of the 
definition of bound and virtual states and show their relation to the parameters introduced 
above. Consider the case when all inelasticities are switched off. Then k e s(E) and K e g(E) 
are given by the simple formulae ()3]) and @. To distinguish between bound- and virtual- 
state situations it is convenient to consider the poles of the denominator E — Ex + Gx{E) 
in the fc-plane, where k = \j2[i(E — E R ). For Ex < E R these poles are 



2/j, p (E r - E x ), k 2 = -id 2[i p (E R - E x ) - ig a bH P - 



(71) 
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The first pole is situated on the first sheet of the genuine two-body Riemann surface, and 
corresponds to the bound state with the energy Ex, and the second pole is at the unphysical 
sheet. The scattering length in the ab channel is positive for this case. On the contrary, for 
Ex > Er the poles are: 



h, 2 = -i 9 -^ ± J^- - 2^{E X - E R ). (72) 

Both poles are located on the unphysical sheet, and the scattering length is negative. For 
dab^p > 8(J5x — Er) the poles are on the imaginary axis, and we deal with a virtual state. 
For small inelasticities, the scattering length acquires a small imaginary part, so that we 
distinguish between the bound and the virtual states by the sign of the real part of a. Thus 
we plot in Figs. |6j|9] the differential branching ratios fl69|) and f JTOj) for both the S- and P- 
wave case and Ex < Er (bound state) and Ex > Er (virtual state). In particular, for the 
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bound-state case we take 

E R = 7MeV, gab = 0.1, E x = E R - 0.5 MeV, T = 1.5 MeV. (73) 



This set of parameters is similar to the one obtained in Ref. 13) from the fitting Belle data 
on the D°D°n and 7t + 7t~J/i/j decay modes of the X(3872) particle. 

As one can see from Figs. |6]and[Hl for an increasing width of the constituent, the lineshape 
in the a[cd] channel starts to develop a peak below the nominal threshold, since the finite 
width of the unstable constituent allows for an overlap with the subthreshold pole. If the 
width is not large, the lineshape above threshold is not affected while, with an increase 
of the width of particle b, the structure present below the nominal ab threshold starts to 
grow very quickly, and it becomes broader. Indeed, already for a width of 1 MeV, the 
structure below the nominal threshold has grown in magnitude and width so much that the 
whole lineshape is just a single broad bump, however, with a noticeable asymmetry letting 
it deviate significantly from the standard Breit-Wigner form. The structure for the inelastic 
channel decreases rapidly with the increase of the width of particle b. The mentioned kind 
of asymmetry, found experimentally in the spectral distribution for the V(4660) observed in 



iP'tttt |8], was interpreted in Ref. 16|], fully in line with the reasoning presented in this work, 
as a clear signal for a dominant molecular /o(980)^' component in the V(4660). In addition, 
the same kind of asymmetry was predicted to show up from an analogous r]' c f (980) bound 



system in the r]' c 7rir invariant mass distribution, yet to be found [17]. 

For the virtual-state case we take the same values for the g a b and r Q as for the bound- 
state case — see Eq. f!73|) — but put E x = E R + 0.5 MeV. Then, as seen from Figs. [7] and 
[9j the effect of the finite width on the a[cd] lineshape is less dramatic since, if the X is a 
virtual state, no pole develops below the threshold. This is already an unambiguous decisive 
difference between the bound state and the virtual state. The inelastic rate decreases with 
the increase of the width of particle b, and the cusplike lineshape, present in the zero-width 
limit, smoothens, so that the shape for the virtual state becomes hardly distinguishable from 
that for the bound state, however, with a shifted peak. 

It follows from the reasoning of Ref. [9[ that the height of the hump above the nominal ab 
threshold contains important information on the nature of the state. We would therefore like 
to repeat here the argument given in Refs. Q, Q that it is dangerous to apply a kinematic 
fit to the data, which moves all relevant signals to above threshold: in this way the signal 
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above the nominal threshold would be risen at the cost of the signal below the threshold. 
Therefore, this procedure might produce data that call for a very different interpretation of 
the X state than that for the original data set. 



VII. INTERFERENCE 



In this chapter we discuss one more effect which plays a role for the lineshapes, namely 



the interferences in the final states 18[. For example, the X(3872) decays in equal parts 
via D°D°* -> D [D°n ] or D°D°* -> D°[D ti } and, clearly, the final states of both decay 
chains interfere. 

If taken naively, the interference can be taken into account through the modification [see 
Eq. d3ZD] 

C(P, Q) -»> C(p, q) = A=(C(p, q) + C(p, q)) (74) 
in Eq. ( 1401) . where p and q are known linear combinations of the momenta p and q: 



m a „ _ m a {m a + m c + m d )m d „ ,._ 

P= P + Qi <? = — <? + 7 p " 75 

m c + m d m c + m d (m c + m d ) z 



The phase-space integration measure is invariant under the variable change {p, q} — > 
{p, q} and the modified Eq. (HDI) reads 

W(M) = j d 3 pd 3 q [|C(p, q)| 2 + Re(C(p, q)C*(p, q))] 5 (M — m a — m c — m d — T 3 (p, q)) , 

(76) 

where the second term in the square brackets describes the interference. It is clear then 
that the interference simply produces an extra contribution to be added to k e g(E), which 
we denote as ki Qt (E) and which can be extracted from Eq. (1761) in the form 

I d 3 pd 3 q Re{C(p, q)C*{p, q))5 (M-m a -m c -m d - T 3 (p, q)) 

k int (E) = k c z{E) J - . (77) 

/ d 3 pd 3 q\C(p, q)\ 2 5 (M — m a — m c — m d — T 3 {p, q)) 

Therefore, in leading order in the interference, the a[cd] rate can be written in the form of 
Eq. (169 p but with the k e ^{E) in the numerator substituted by the sum k e g(E) + k int (E): 

dBi(a[cd]) = B _L g ab (k eS (E) + k int (E)) 

dE ~ 27r(E-E x + Sfn eS (E)y+ 1 j(r + g ab k cS (E)f 
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FIG. 10: The a[cd] state production rate for the S-wave resonance for bound-state case (the first 
plot) and for the virtual-state case (the second plot). Solid line: no interference; dashed line: naive 
interference; and dotted line: improved interference. As before, we set B = 1. 

Clearly, taking the interference into account in the numerator only while neglecting it in 
the denominator (higher orders in the interference) leads to a violation of unitarity, which 
can be cured by solving the full coupled-channel problem with the effects of dynamical 
building blocks of the constituents taken into account from the very beginning. However, 
bearing in mind a very moderate effect of the interference over the lineshapes, we resort to a 
much simpler approach which amounts to taking ki nt (E) into account both in the numerator 
and in the denominator simultaneously: 



We call this an improved interference. 

In Fig. [TD]we plot the rate dBr(a[cd])/dE for an S"-wave resonance and for the following 
three cases: (i) no interference [as given by Eq. ( 1691) ]. (ii) naive interference [as given by 
Eq. (1751) ]. and (hi) improved interference [as given by Eq. (1791) ]. For definiteness, we fix all 
parameters compatible with the X(3872) charmonium, so that the parameters used are 

E R = 7 MeV, g ab = 0.1, E x = E r t 0.5 MeV, T = 1.5 MeV, r R = 100 keV, (80) 

where upper (lower) sign in Ex corresponds to bound(virtual)-state case. Finally, we set 
= m c — m o and irtd = m n , with m and m n being the masses of the D° and 7r° mesons, 
respectively. 



dBr(a[cd]) 
dE 




2tt (E - E x + ?fK eS (E)y + \ (r + g ab (k efi (E) + k int (E))) 



1 g ab (k eS (E) + k int (E)) 



(79) 
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The corresponding plots for the P-wave exhibit exactly the same behaviour as for the 
S-w&ve and are omitted here. 

As one could anticipate, the net effect of the interference is an increase of the rate in the 
region of the bound-state peak, while the above-threshold rate remains practically the same. 

A proper account for the interference obviously requires a more rigorous treatment of 
the coupled-channel problem. Indeed, in the case of the X(3872), a potentially important 
piece of interaction is omitted in the Hamiltonian (|16p . which is responsible for the D°D *- 
D°D°* transition, namely the pion exchange. This interaction is often considered as the 
main binding mechanism in the D°D°* system 19l-l21| , and the X is generated dynamically, 
similarly to the deuteron. It was noticed, however, that, as the nominal threshold for D°D°* 
is only about 7 MeV larger than the three-body D°D 7r° threshold, the pion may go on-shell 
and, as a result, such a binding may not be strong enough (see, for example, Ref. |22|). 
It was argued in Ref. 23j that the effect of the ir° exchange for the D°D°* system can 
be treated perturbatively, in contrast to the two-nucleon system. If it is indeed the case, 
the expression ( 1TH1) for the D°D ir rate should be appropriate. This statement is further 
supported by the observation that the rates described with the "naive" and the "improved" 
prescriptions for the interference do not differ a lot. 



VIII. SUMMARY 



In this paper we investigated the effect of the finite width of a constituent of a composite 
state on its lineshape. In a coupled-channel formalism we derived exact formulae for such 
lineshapes and compared them to simple prescriptions found in the literature. We found that 
such simple prescriptions work only for extremely narrow constituents located far away from 
the three-body threshold characterised by A = Pr/(2Pr) <C 1, where Pr denotes the width 
of the constituent and Er its resonance energy, measured with respect to the three-body 
threshold. In this regime the resonance width can safely be treated as energy independent. 
The famous charmonium state A (3872) gives an example of such a safe situation with A 1 
in view of an extremely small width of the D* meson. However, for an increasing width of 
the constituent (and correspondingly an increasing value of A), the exact formulae start to 
deviate significantly from the simplified ansatz. The situation cannot be improved by using 
an energy-dependent width for the constituent instead of a constant width — the three-body 
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threshold behaviour gets spoilt in this case, while the near-two-body threshold behaviour is 
not improved. We argue therefore that the exact formulae have to be employed in such a 
situation. 

The effect of the finite width of a constituent on the lineshape in the elastic channel of the 
composite state X is quite moderate, if the X is a virtual state. However, for the X being a 
bound state, a nonnegligible bound-state peak starts to develop below the nominal threshold 
for an increasing constituent width, as already stressed in Ref. [3[. For sufficiently large 
values of the width, the bound-state peak becomes broad enough to distort the lineshape 
both below and above the threshold, so that the resulting lineshape takes the form of a broad 
hump that one can hardly distinguish in shape from the virtual-state situations. However, 
the peak position appears to be shifted, since for the virtual state the distribution is peaked 
above the nominal two-body threshold while for the bound state it is peaked below. Clearly, 
to distinguish the two experimentally calls for a quite high resolution. 

The inelastic production rate is less affected by the finite-width effects. However, only 
for extremely small values of the width one is able to distinguish between the threshold 
cusp scenario, in case of the virtual state, and a smooth below-threshold peak, in case of 
the bound state. For an increasing width, a smeared broad structure is developed in both 
cases. This makes it difficult to distinguish between the two scenarios, although the peaks 
of the two distributions appear to be shifted also in this case. It is important to stress that 
the distributions turn out to be very asymmetric — inconsistent with a Breit-Wigner shape 
- in both cases. This asymmetry can be taken as a strong signal towards a composite 
two-hadron structure in distinction to a genuine, compact quark state. 

As a related problem, we studied the interference in the final state. For the sake of 
definiteness, we focused on the case of the X(3872) charmonium and fixed all free parameters 
accordingly. The conclusion we find is that, while the effect of the interference on the elastic 
lineshape is indeed sizable for the bound state and below the nominal threshold, in line 
with the findings reported in the literature, this effect is negligible for the bound state 
above the threshold as well as for the virtual state both below and above the threshold. 
However, for an increasing width one is to include the effects of interference. Furthermore, 
the naive inclusion of the interference meets severe problems with unitarity, so that, while for 
intermediate values of the width, an improved interference suggested in this paper should be 
used, for a width large enough one is finally forced to solve the full coupled-channel problem 
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with the effects of dynamical building blocks of the constituents taken into account from the 
very beginning. 

We treat our findings as an important step in building a universal lineshape parameter- 
isation scheme suitable for fitting experimental data for various near-threshold resonances 
and thus we believe that our insights will prove useful for the analysis of future experiments. 
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